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The resummation of soft gluon exchange for QCD hard scattering requires a matrix of anomalous 
dimensions. We compute this matrix directly for arbitrary 2 — » n massless processes for the first 
time at two loops. Using color generator notation, we show that it is proportional to the one- 
loop matrix. This result reproduces all pole terms in dimensional regularization of the explicit 
calculations of massless 2^2 amplitudes in the literature, and it predicts all poles at next-to- 
next-to-leading order in any 2 — > n process that has been computed at next-to-leading order. The 
proportionality of the one- and two-loop matrices makes possible the resummation in closed form 
of the next-to-next-to-leading logarithms and poles in dimensional regularization for the 2 — > n 
processes. 



The calculation of high-energy cross sections in per- 
turbative quantum chromodynamics (QCD) for hadronic 
collisions involves the factorization of long- and short- 
distance effects. Sensitivity to long-distance dynamics 
is enhanced by powers of logarithms whenever there is 
an incomplete cancellation of parton emission and vir- 
tual corrections. In such situations, it is useful to or- 
ganize, or resum, these corrections to all orders in per- 
turbation theory. Correspondingly, in partonic scatter- 
ing or production amplitudes, it is necessary to organize 
poles in e that arise in dimensional regularization (with 
D = 4 — 2e) . The resummation of these poles and related 
logarithmic enhancements is well-understood for inclu- 
sive reactions mediated by electroweak interactions, such 
as the Sudakov form factor 0, 0] and in Drell-Yan pro- 
cesses yj- With recent advances in the computation of 
splitting functions many such corrections can be re- 
summed explicitly to next-to-next-to-leading level. Their 
structure at arbitrary level is known to be determined by 
a handful of anomalous dimensions. 

The situation for QCD hard scattering processes con- 
taining four or more partons — critical to understanding 
many types of backgrounds to new physics at hadron 
colliders [{| — is more complex. Resummation beyond 
leading logarithms or poles requires a matrix of addi- 
tional anomalous dimensions [y, 151 la 13 • These matrices 
are found in turn from the renormalization of the vacuum 
matrix elements of products of Wilson lines, one for each 
external parton in the underlying process |jj ■ In this pa- 
per, we investigate the structure of the two-loop anoma- 
lous dimension matrix. We will find that, remarkably, 
for every hard-scattering process involving only massless 
partons, this matrix is proportional to the one-loop ma- 
trix. We will concentrate below on the role that the ma- 
trix plays in partonic amplitudes. The full calculation of 
the two-loop matrix will be given elsewhere El ■ ^ n this 
paper, we provide the simple calculation that is at the 
heart of the main result. We will show that certain color 
correlations due to two-loop diagrams that couple three 
Wilson lines vanish identically. We will also provide an 
explicit expression in terms of color generators [TTL IT^ 
for all single-pole terms in massless 2 — > n amplitudes. 



We consider a general process involving the scattering 
of massless partons, which we denote by "f ' : 



f: fi{pi,ri) + f 2 {P2,r 2 ) 

-> fs(p3, ra) H h / n +2(Pn+2,r n+2 ) ■ 



(1) 



The fi are the flavors of the participating partons, which 
carry momenta {pi} and color {r^}. Adopting the color- 
state notation of Ref. we represent the amplitude 
for this process as \M{). 

It is convenient to express these amplitudes as vectors 
with C elements in the space of color tensors, for some 
choice of basis tensors {(c/)r i} 0, El Eli 
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We will analyze these amplitudes at fixed momenta pi 
for the participating partons, which we represent as 
Pi = QPi , Pi = , where the Pi are four- velocities, 
and where Q is an overall momentum scale. 

In dimensional regularization, on-shell amplitudes may 
be factorized into jet, soft and hard functions that de- 
scribe the dynamics of partons collinear with the exter- 
nal lines, soft exchanges between those partons, and the 
short-distance scattering process, respectively. This fac- 
torization follows from the general space-time structure 
of long-distance contributions to elastic processes [3, [w| . 
The general form of the factorized amplitude , fo r equal 
factorization and renormalization scales fi, is 13] 
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where JW is the jet function for external parton i, Sf 
is the soft function, and H{ is the hard (short-distance) 
function. 
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The jet function for parton i can be expressed to all 
orders in terms of three anomalous dimensions, K\ l \ 
and Yk i °f which the first is determined order-by-order 



from the third. The general form of the jet function, and 
its expansion to second order is given by (expanding any 
function as f(a s ) = E„ f (n) ) 0, * 
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In the expansion we use the Z3-dimcnsional running-coupling, evaluated at one-loop order, 



a s (/i,e) = a s (» f ) 
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with the one-loop coefficient /3q = 11CU/3 — AT F n F /3. The corresponding expression for the soft matrix is 



S f ( — ^ 3 ,a s (fi),e ) = P exp 
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, (6) 



where uq — fi 2 /Q 2 , so that pi ■ Pj/uo = Sij/^i 2 . 

Expanding Q W = + e{/M' + . . . , one finds from 
Eq. J3J| the single pole in e in the logarithm of the jet 
function at two loops. For the quark case this term is 
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using values of Q^(e) from ref. ^^|. Notice the contri- 
bution from the running of the finite term at one loop, 
which appears as an 0{e) contribution in (jMf 1 ). 

The one-loop soft anomalous dimension in color- 
generator form is 



111 



(8) 



where Sjj = (pi + pj) 2 , with all momenta defined to flow 
into (or out of) the amplitude. The Tj are given explic- 
itly by color generators in the representation of parton 
i, multiplied by ±1: plus one for an outgoing quark or 
gluon, or incoming antiquark; minus one for an incoming 
quark or gluon, or outgoing antiquark. The color gen- 
erator form for the anomalous dimension matrix is more 
flexible, but less explicit, than the corresponding matrix 



expressions in a chosen basis of color tensors for the am- 
plitude. An example of the latter for qq — ► qq scattering, 
in an s-channel t-channel singlet basis, is 



• s f 



i 



(U - T) + 2 C F S 
(T-U) 



(S-U) 
(U - S) + 2 C F T 



where T = In 



(9) 

and so on for the other Mandclstam 



invariants, defined by s — S12, t — S13, u — S14. Re- 
summed cross sections are determined by the eigenvalues 
and eigenvectors of these matrices 0, . 

We are now ready to provide our result for the full 
two-loop soft anomalous dimension matrix, 
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K 
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Here K = CU (67/18 - C(2)) - lOT F n F /9 is the same con- 
stant appearing in the relation between the one- and two- 
loop Sudakov, or "cusp" anomalous dimensions 0, : 
7 M = 2Cj(a s /7r)[l + (a a /w)K/2]. Remarkably, relation- 
ship I|1U|) holds for an arbitrary 2 — > n process, even 
though the two- loop diagrams shown in Fig. ^ apparently 
couple together the color factors of three eikonal (Wilson) 
lines coherently. We derive Eq. IjlUI) using the color gen- 
erator formalism; however, the result is completely gen- 
eral, and applies to explicit matrix representations such 
as Eq. 10. 

Following the method described in detail at one loop 
in Ref. Q, and extended to two loops in Ref. t ne 
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FIG. 1: Two-loop diagrams involving three eikonal lines, 
two-loop anomalous dimension is found from the residue 



of single-pole terms in suitable combinations of Wilson 
lines computed at two loops. The simplicity of the re- 
sult 1)10(1 follows from the special properties of the di- 
agrams of Fig. ^ which connect three different eikonal 
lines. First consider Fig. QJl, where a three- gluon cou- 
pling ties together three eikonals labelled va, vb and vc, 
in an otherwise arbitrary eikonal process. We shall prove 
that this integral is zero, as long as we take the eikonals 
va and vb to be lightlike. 

Since va and vb are lightlike, we choose a frame in 



which v^ = 8^ and 



<5 M+ . Components transverse 



to va and vb carry the subscript T. The eikonal integral 
in momentum space is then 



Fia(vA,VB,vc) = I d D kid D k 2 



k\ + ie k\ + ie (ki + k 2 ) 2 + ie k x + ie k 2 + ie 
v c (k+ - k+) + v+ (fcf - k 2 ) - v c ,t ■ (ki, T ~ k 2 , T ) + «c { k i + 2k 2 ) + v c (~ 2k t - k t) 



u c (fcf + k%) + Vq (fcj + k 2 ) - v c ,t ■ (ki,T + k 2 ,r) + ■ 



(11) 



where the term in square brackets is the three-gluon ver- 
tex momentum factor. We now introduce a change of 
variables (with unit Jacobean) from momenta kf to kf, 



(k^,k 1 ,k 1 , T ) = ((k 2 ,( 1 k 2 h ,k 2 , T ), 
(fe^,fcj,fc 2 ,r) = (CK , C" 



(12) 



where C = v c l v c • ^ provides an expression identical to 
Fia{vA, vb, vc), Eq. ifTTJl. but of the opposite sign. The 
integral corresponding to Fig. ^ therefore vanishes. 

Regarding Fig.^ 3 , the same change of variables yields 
l/[(v c ■ (h + k 2 )){v c ■ kx)] = l/[(v c ■ (h + h)){vc ■ h)}, 
from which it is easy to show that this diagram reduces 
to the product of one-loop diagrams, and so does not 

contribute to the two-loop anomalous dimension. In- 

( 2 ) 

deed, the only nontrivial contributions to at two 
loops involve only two eikonal lines. Using results from 
refs. [TtI ITsj. the color structure of these contributions 
reduces to that of a single gluon exchange. The sum of 
the diagrams then modifies the one-loop result by the 
same multiplicative factor as for the cusp anomalous di- 
mension, which gives Eq. (|1(J|) . 

The explicit expression for single poles in 2 — > n am- 
plitudes is easily found from Eqs. Q and © using the 
explicit form of the two- loop matrix (|10f) . 
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Here we normalize the one-loop hard scattering by ab- 
sorbing into it all finite terms from the jet functions, or- 
der by order, and \H^\q)), \Ti!^ '{{))) are this function 
and its derivative with respect to e, respectively, evalu- 
ated at e = 0. (This absorption is possible to any loop 
order because the jets are diagonal in color.) Explicit 
comparison 0] shows that this simple result agrees with 
all single-pole terms found at 2 — > 2 in the literature, as 
summarized for example in Refs. [l9ll20| . It also predicts 
all such poles in a 2 — > n process, once the one-loop hard 



part 



(i) 



is known. 



We remark that an analogous anomalous dimension 
matrix for Wilson lines has been computed in Ref. [l8| . 
in the forward limit t — ► 0. This limit is a singular one, 
with respect to our arguments regarding Fig.^ thus our 
results and theirs are not directly comparable. 

It is also worth remarking on the relationship between 
our results and the influential alternative formalism of 
Ref. [T^ . in which both pole and finite terms are put 
into an exponential form to two loops. We may think 
of these as alternative schemes for organizing infrared 
poles. When explicit calculations are organized according 
to the scheme of Ref. ^3 > more complex color products 
are found, namely i^T^T^Tl = -[T* • Tj , T 3 • T fc ], 

appearing in the matrix H^ 2 ' at order 1/e 0,[2(|. Such 
products are not encountered in the resummation scheme 
described above. These differences in color structure, 
however, are by no means disagreements. They arise 
from a particular commutator, between the one-loop fi- 
nite terms that are exponentiated in the formalism of 
Ref. ^3 , an( i the one-loop soft anomalous dimension 
(i) 



matrix T K g' . The result of performing the commuta- 
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tor [lj} agrees with the form of H^ 2 ) in Ref. for 
2 — ► 2 processes, and with that proposed in Ref. [2l| 
for 2 gluon — > n gluon processes, based on consistency 
of collinear limits. 

Given an explicit two-loop amplitude, the strategy de- 

(2) 

scribed here may be reversed, and Tg extracted di- 
rectl y fr om the amplitude. This approach was adopted in 
Ref. |22| for the case of quark-quark elastic scattering, in 
the context of electroweak Sudakov corrections. The orig- 
inal version of Ref. [12 differs from ours due to omission 
of the commutator contribution described above. The 
authors have informed us that a revision is in prepara- 
tion. 

Similar remarks apply to the color-diagonal single 
poles given in Eq. These coefficients do not equal 

(2) 

the corresponding coefficients in the formalism of 

Ref. ^3 ' but they are connected |23j | . The difference can 
be related precisely to the different treatment of finite 
terms in the two approaches fTo| . 

In addition to clarifying the structure of singular terms 
in calculations of 2 — > 2 processes at two loops, the re- 
sults outlined here have potentially useful consequences 
and suggest further directions of research. Eq. i|13|) pre- 
dicts the two-loop pole structure for any 2 — > n process, 
in color-generator form, for any process whose one-loop 
hard function is known to 0{e). 

Another practical consequence is that, because the 
one- and two-loop anomalous dimensions are propor- 
tional, all terms in the expansion of the soft function 
commute to next-to- next-to-leading level (NNLL), and 
in this approximation, the ordering operator P can be 
dropped in Eq. (JSJ). Thus, once the color eigenstates of 



the one-loop matrix are known, the same states will diag- 
onalize the two-loop matrix. A semi-numerical approach, 
bypassing diagonalization, is to simply exponentiate the 
relevant matrices in any convenient basis [9j. Given the 
relation H10|) . this is now possible at NNLL as well as 
NLL. 

The study of these matrices for processes beyond 2 — > 
2, already begun in Ref. [24J, is clearly an important chal- 
lenge. Another intriguing question is whether the propor- 
tionality HI Up might extend beyond two loops, whether in 
QCD or any of its allied gauge theories. If so, it could 
have consequences for the interpretation of infrared di- 
verences in the relevant theory. The extension, and/or 
modification of the results above for the production of 
massive colored particles is another important direction 
for research. 
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